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Every  nonsingular  totally  positive  m-banded  matrix  is  shown  to  be  the 
product  of  m  totally  positive  one-banded  matrices  and,  therefore,  the  limit 
of  strictly  m-banded  matrices.  This  result  is  then  extended  to  (bi) infinite, 
'nonsingular' ,  totally  positive  matrices.  In  the  process,  such  matrices  are 
shown  to  possess  at  least  one  diagonal  whose  principal  sections  are  all 
nonzero. 
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SIGNIFICANCE  AND  EXPLANATION 


This  report  deals  with  certain  technical  questions  regarding  totally 
positive  matrices,  i.e.,  matrices  with  all  minors  nonnegative.  Such  matrices 
make  their  appearance  in  various  areas  of  analysis.  Particular  emphasis  here 
is  on  totally  positive  matrices  which  are  also  banded  (meaning  that  all 
nonzero  entries  occur  in  just  a  few  consecutive  bands  or  diagonals)  and 
biinfinite. 

Such  matrices  occur,  for  example,  in  the  study  of  standard  spline 
approximation  schemes  such  as  spline  interpolation  at  knots  or  least-squares 
approximation  by  splines.  Certain  questions  about  such  a  matrix  are  much  more 
easily  answered  when  it  is  known  that  the  matrix  is  strictly  banded.  This 
means  that  all  the  elements  in  the  two  outermost  nontrivial  bands  of  the 
matrix  are  nonzero.  The  answers  can  then  also  be  applied  to  a  totally 
positive  matrix  which  is  merely  banded  provided  one  knows  that  there  are 
strictly  banded  totally  positive  matrices  arbitrarily  close  to  the  matrix  of 
interest.  This  useful  fact  is  proved  in  this  report  by  showing  that  a  non¬ 
singular  totally  positive  matrix  can  be  factored  into  certain  bidiagonal 
totally  positive  matrices,  even  when  the  matrices  involved  are  biinfinite. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sunsary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A  FACTORIZATION  OF  TOTALLY  POSITIVE  BAND  MATRICES 


Carl  de  Boor  and  Allan  Pinkus 

1.  Introduction*  In  this  paper,  we  prove  the  result,  needed  in  [1J,  that  a  totally 
positive  biinfinite  band  matrix  is  the  limit  of  ' strictly  banded’  totally  positive  matrices 
of  the  same  band  type.  But  the  tool  developed  for  the  proof,  viz.,  the  factorization  of 
such  matrices  into  'one-banded*  totally  positive  matrices,  is  of  independent  and,  perhaps, 
greater  interest. 

We  came  to  consider  such  factorizations  because  of  the  recent  paper  by  Cavaretta, 
Dahnten,  Micchelli  and  Smith  [3]  in  which  such  a  factorization  is  derived  for  strictly 
banded  totally  positive  matrices.  But  we  were  unable  to  adapt  their  arguments,  which 
involve  limits  of  ratios  of  entries  in  a  certain  matrix  inverse,  to  our  situation,  and 
ended  up  constructing  the  needed  factors  by  the  more  familiar  device  of  elimination 
instead.  The  factorization  is  first  established  for  finite  matrices  and  is  then  extended  to 
biinfinite  matrices  by  a  limiting  argument.  For  this,  we  found  it  necessary  to  first  prove 
that  such  a  totally  positive  'nonsingular*  matrix  has  at  least  one  diagonal  with  the 
property  that  all  square  finite  sections  which  are  principal  for  that  diagonal  are 
nonsingular. 

2.  Bandedness.  The  r-th  diagonal  or  band  of  a  matrix  A  is,  by  definition,  the 
sequence  (A(i,i+r))  .  As  in  [2],  we  call  a  matrix  A  m- banded  if  all  nonzero  entries  of 
A  can  be  found  in  at  most  m+1  consecutive  bands.  Explicitly,  the  matrix  A  is  m-banded 
if 

for  some  i  ,  A(i+l,j)  ^  0  implies  i  <  j  <  i+m  . 

If  both  1  and  m-1  are  nonnegative,  then  the  m+1  nontrivial  bands  include  the  ‘main 
diagonal'  or  zeroth  band,  with  1  bands  to  the  left  of  it  and  r  :«  m-i  bands  to  the 
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right  of  it.  In  this  situation,  we  will  at  times  call  such  a  matrix  more  explicitly 
(4,r)-banded. 

We  call  a  band  matrix  strictly  banded  if  the  leftmost  and  the  rightmost  nontrivial 
band  contains  no  zero  entries.  Among  banded  matrices,  the  strictly  banded  ones  are 
particularly  easy  to  treat  since  they  correspond  to  nondeqenerate  difference  operators. 


3.  Total  positivity.  A  matrix  is  said  to  be  totally  positive  (or  TP)  in  case  all  its 
minors  are  nonnegative. 

We  use  the  abbreviation 

:m  (A(wC  v-i 

for  the  sxt  matrix  which  has  its  (y,v)-entry  equal  to  A(iy,jv)  .  Further,  if  I  and  J 


are  index  sets,  then 


i.,  •••  ,1 


4,  *tj]  XL’  /] 

IfJ  J  fJ 


with  i1#  ...,  ig,  and  jy  ...,  jt  the  elements  of  I  and  J  ,  respectively,  in 
increasing  order.  Occasionally,  we  will  use  the  additional  abbreviation 

a(iJ  !«  a[*]  . 

Finally,  replacing  the  square  brackets  by  round  brackets  gets  us  from  the  matrix  to  its 
determinant: 

A(  •”  )  :*=  det  A[  ] 

•  •  •  »  *  • 

We  will  make  repeated  use  of 


Sylvester's  determinant  identity  (SOI).  If.  a(*]  i  0  and  B  is  the  matrix  obtained 


from  A  by 


then 


B(  i,  j)  A(|JJ  .  all  (i,j)  t  I*J  , 

■(j:j  -  • 


The  submatrix  Aj  j  is  called  the  pivot  block  since  the  identity  is  proved  by 


observing  that  B  is  the  Schur  complement  of  A^j 


,  i.e.,  the  interesting  part  of  what  is 


left  in  rows  \I  and  columns  \J  after  rows  I  have  been  used  to  eliminate  variables 
J  from  the  other  rowsi  see,  e.g.,  Gantmacher  [5;p.31]  or  Karlin  [ 7 ; p. 3 ]  .  In  particular. 

Corollary  1  •  rank  a[  * , U  *  ]  “  1 1 1  +  rank  B  l.V.1  HSL  X-nI-0-J-nJfc 

J  U  J  J 

Corollary  2.  B  is  again  TP  if  A  is. 

Another  result  which  may  be  proven  by  Sylvester's  determinant  identity  (using 
induction  on  |l|»  see  Gantmacher  <.  Krein  [6ip.108]  or  Karlin  [7;p.88]  )  is 

Hadaaard's  inequality.  If  A  is  TP  and  I  *  I*  u  I"  with  I'n  I"  =*  0  ,  then 

A( I )  <  A(I')A(I")  . 


4.  Shadows.  In  this  section,  we  prove  an  ancillary  result  concerning  the  existence  of 
a  diagonal  in  a  TP  biinfinite  matrix  which  could  serve  as  the  main  diagonal  in  a  triangular 
factorization,  i.e.,  a  diagonal  all  of  whose  principal  sections  are  nonsingular. 

A  zero  entry  in  a  TP  matrix  usually  'throws  a  shadow*.  By  this  we  mean  that  usually 
all  entries  to  the  left  and  below  it,  or  else  all  entries  to  the  right  and  above  it,  are 


ino 
^0 

A(ip,Jo)  and,  correspondingly,  call  the  submatrix 
A(i(j,jQ).  Then  the  following  lemma  is  known. 


also  zero.  More  precisely,  call  the  submatrix  a[ 


*1 


the  left  shadow  of  the  entry 
i«i 

]  the  right  shadow  of 

3  30 


mu  If  A  is  TP  and  A(i0,j0)  -  0  ,  but  neither  A(*,j0)  nor  A(  Iq  ,  - )  is  zero, 
then  either  the  left  or  the  right  shadow  of  A(ig,j0)  is  zero. 


Proof.  By  assumption,  A(ig,jj)  ¥  0  for  some  .If  <  jg  ,  then  the  right 

shadow  of  A(iQ,j0)  can  be  seen  to  be  zero  as  follows.  First,  for  any  i  <  iQ  , 

i  rin 

0  <  A(  U)  -  -  A(i  ,J  )A(i,j  >  <  0 

3,»30  0  1  u 
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and  Adg.jf)  f  0  Implies  that  A(i,jQ)  »  0  for  all  i  <  1q  .  Hence  there  then  exists  i., 


>  i0  for  which  Ad^.jg)  /  0  .  But  now,  for  any  i  <  i  ,  j  >  j  , 

J  i  ,i  00 

0  *  A^j  “  ‘  A<  WA(i,j)  4  0 

and  Adj/jg)  /  0  implies  that  A(i,j)  -  0  . 

Finally,  if  instead  >  jQ  ,  then  the  left  shadow  of  A(ig,jQ)  is  similarly  seen  to 
be  zero.  |  |  | 

As  an  application  for  later  use,  note  that  a  zero  in  the  lower  triangular  part  of  an 
invertible  TP  matrix  necessarily  throws  a  left  shadow  since  all  diagonal  entries  are 
nonzero,  by  Hadamard's  inequality. 

More  generally,  for  any  section  of  A  ,  i.e.,  any  submatrix  Aj  ^  of  A  made  up  of 

consecutive  rows  and  columns  of  A  ,  we  call  the  submatrix  of  A  having  hj  j  as  its 

upper  right  corner  the  left  shadow  of  >  and,  correspondingly,  we  call  the  submatrix 

having  *I#J  as  its  lower  left  corner  the  right  shadow  of  Aj  j  .  Then  we  have  the 
following  generalization  of  the  lemma. 

Proposition.  Xt_  a  is  TP  and  Aj  j  is  a  singular  section  of  order  n  and  rant 

n-1  ,  while  both  a[*]  and  A^]  are  of  full  rank  n  ,  then  either  the  left  or  the  right 

shadow  of  A.  .  has  rank  n-1  . 

-  x ,  j  - 

Remark.  As  in  the  case  n-1  discussed  earlier,  we  will  describe  this  last  situation 
by  saying  that  such  a  section  Aj  j  'throws  a  (left  or  right)  shadow'. 

Proof.  By  assumption,  we  can  choose  (ig,jg)  £  IXJ  80  that  Alj.)  *  0  >  with  I*  :  = 

I  \  {iQ}  ,  J’  J\{jg}  .  The  assumptions  imply  that  the  Schur  complement  of  Aj,  j,  ,  i.e., 
the  matrix  B  given  by 

B(r,s)  A([aj  y )  *  a11  d,s>  4  I'*J'  , 

is  again  TP  (by  Corollary  2  of  SDI)  and  vanishes  at  (i0,j0)  while  (by  corollary  1  of  SDI) 
neither  B(iQ,*)  nor  B(*,jg)  is  zero.  The  lemma  therefore  implies  that  either  the  left 


k 


1 


or  the  right  shadow  of  B(ig/jQ)  is  zero,  and  Corollary  1  of  SDI  then  implies  that  either 
the  left  or  the  right  shadow  of  Aj  j  has  rank  n-1  .  | | | 

IfXN 

Corollary,  _I£  A  is  an  infinite  TP  matrix,  e.g»,  A  t  R  ,  then  all  rows  and  all 
columns  of  A  are  linearly  independent  if  and  only  if  A(I)  ?  0  for  all  finite  I  N. 

Proof.  If  A(I)  =0  for  some  finite  I  £  N,  then  Hadamard' s  Inequality  implies  the 
existence  of  some  n  e  N  for  which  A(1,...,n)  =  0  while  A(1,...,n-1)  f  0.  The 
proposition  then  implies  that  either  the  first  n  rows  or  else  the  first  n  columns  of 
A  are  linearly  dependent.  | | I 

We  now  state  and  prove  the  corresponding  result  for  a  biinfinite  TP  matrix.  This  is 
somewhat  harder  since  it  is  not  clear  a  priori  which  band  is  to  play  the  role  of  main 
diagonal . 

We  concentrate  on  principal  sections  for  a  band:  A  principal  section  for  band  r  is 
any  submatrix  of  the  form  Aj  I+Jf  ,  with  I  an  interval.  In  other  words,  such  a  principal 
section  for  band  r  is  (i)  square,  (ii)  made  up  of  consecutive  rows  and  columns,  and  (iii) 
has  a  piece  of  band  r  as  its  main  diagonal.  We  call  such  a  principal  section  minimally 
singular  if  it  is  singular  but  contains  no  smaller  principal  section  for  the  same  band 
which  is  also  singular.  Note  that,  by  Hadamard' s  inequality,  every  principal  section 
containing  a  singular  one  for  the  same  band  is  itself  singular. 

Theorem  A.  Let  A  be  biinflnlte  TP  and  assume  that  not  all  minimally  singular 
sections  of  A  throw  their  shadow  in  the  same  direction.  Then,  all  rows  and  all  columns 
of  A  are  linearly  Independent  if  and  only  if  all  principal  sections  for  some  band  are 
nonsingular,  l.e., 

there  exists  r  so  that,  for  all  intervals  I  ,  a(  t  0  . 
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*> 


» 


Proof.  The  sufficiency  of  the  condition  is  obvious.  So,  assume  that  all  rows  and  all 
columns  of  A  are  linearly  independent.  If  every  band  r  has  a  singular  principal 
section,  then  every  band  has  a  minimally  singular  one  and,  all  rows  and  all  columns  being 
linearly  independent,  each  of  these  throws  a  shadow,  by  the  proposition.  By  assumption,  not 
all  of  these  shadows  go  in  the  same  direction,  i.e.,  there  exist  bands  r  and  s  so  that 

some  principal  section  for  r  throws  a  shadow  to  the  left  while  some  principal  section 

for  s  throws  a  shadow  to  the  right. 

We  may  assume  that  r  <  s  .  For,  if  the  left  shadow  is  of  rank  L  ,  say,  then  every 

band  q  <  r  has  a  minimally  singular  principal  section  (of  order  <  L+1 )  inside  that 

shadow,  and  each  must  throw  its  shadow  to  the  left  since  otherwise  the  union  of  the  two 

shadows  would  contain  a  strip  of  width  >  L+1  and  of  rank  <  L,  thus  contradicting  the 

linear  independence  of  rows  (or  columns)  of  A  . 

Further,  since  every  band  has  a  shadow  throwing  section,  we  may  assume  that  s  »  r+1  . 

More  explicitly  now,  we  assume  that,  for  some  r  ,  there  is  L  >  0  and  i  so  that 

A  s-  a[  *  '  is  minimally  singular  with  its  left  shadow  of  rank  L  while, 

r  i+ 1  §  •  •  •  t  i+  r+ l 

r  fc-R  see  i 

for  some  R  >  0  and  some  k  ,  A  . ,  :=  a[  '  ’**  is  minimally  singular  with 

^  *  K'  r+i •••  jK+r+i 

its  right  shadow  of  rank  R  . 

There  are  three  cases,  depending  on  the  relative  position  of  these  two  submatrices. 

r  i  •  K  i 

Case  1:  k-i  >  max{L,R}  .  Then  all  columns  j  in  A[  '  ***  '  J  with  j  <  i+r  are 

i+r  i+r+L 

.  +  + 

i  ♦  x  x  x  x  x 

x 

x 

x 

x 

xxxxxx+k 
k+r+1-R  k+r+1 

linear  combinations  of  the  L  columns  i+r+),  i+r+L,  while  all  columns  j  with 

j  >  k+r  are  linear  combinations  of  the  R  columns  k+r+1-R,  ...,  k+r  ,  and  this  holds 
even  if  ,  e.g.,  L  ■  0  since  then  all  columns  j  with  j  <  i+r  are  zero  themselves. 

Since  i  +  r+1  <  k  +  r+1  -  R  and  i+r+L<k+r,  all  columns  are  linear  combinations 
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**•  • 


of  the  k-i  columns  i+r+1,  . k+r  .  We  conclude  that  a[*'  ,lc]  has  only  rank 
k-i,  a  contradiction  to  the  assumed  linear  independence  of  all  rows* 

Case  2s  i+r+L  -  (k+r+1-R)  >  max{L,R}  .  This  is  treated  analogously.  It  leads  to 
columns  k+r+1-R,  . ..,  i+r+L  being  dependent*  again  a  contradiction. 

This  leaves 

Case  3:  max{k-i,  i+L  -  (k+1-R)}  <  max{L,R}  .  In  this  case  i-k  <  1  +  max{L,R) 

(L+R)  “  1  -  min{L,R}  <  and  so 

min{L,R}  <  k-i  <  max{L,R}  . 

We  claim  that  this  contradicts  the  minimality  of  the  two  sections  Ar  and  Ar+1  chosen. 
Assume  without  loss  that  L  <  R  .  Then,  since  L  <  k-i  <  R  ,  the  r-band  section  lies 

inside  the  (larger)  (r+1)-band  section  Ar+^. 


i+r  i+r+L 


i  ♦ 
i+L  + 


X  !*x*  "  X I  • 

1  l 

x  l*.*-.!  * 

XXX* 


k-R 


SDI 


k+r+1-R 


k+r+1 


©  ®  •  • 

©  ©  •  • 


By  the  minimality  of  A  ..  ,  the  L-section  a[  *  '  ***  ,i*1'  )]  principal  for  band 
r+1  is  nonsingular.  In  its  Schur  complement  B  ,  the  section  Ar  appears  as  a  zero  which, 
by  assumption,  throws  a  left  shadow.  But  this  zero  appears  on  the  next-to-main  diagonal  of 
the  submatrix  Br+^  of  B  corresponding  to  Ar+1  and  this  implies  det  Br+1  to  be  the 
product  of  two  of  its  proper  complementary  minors.  But  then,  by  the  singularity  of  Ar+1  , 
hence  of  Br+1  ,  some  proper  principal  section  of  Ar+^  must  be  singular,  contradicting 
the  minimality  of  Ar+1  .  | | ( 
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The  assumption  that  there  ought  to  be  shadows  thrown  in  both  directions  cannot  be 
ommitted,  as  the  following  example  shows.  Let 


/ 

'  1/<i-j)  1  , 

B(i,  j)  :=  < 

0 

i<3 

* 

Then  B  is  TP  and,  with  i1  <  ...  <  is. 

V 

• 

V 

V 

•r-» 

' 

we  have 

^ /  •••  i 

*L'  *)  > 

. 3s 

0  iff  i 

>  it 

,  all  t  • 

From  B  ,  construct  A  by  deleting  every  other  row.  More  explicitly, 

A( i, j)  i=  B( 2i ,  j )  ,  all  i,j  . 

Then  A  is  again  TP  and  biinfinite  and  its  rows  and  columns  are  linearly  independent.  But 
now  every  band  r  has  singular  principal  sections;  for  example,  A(i,i+r)  =  0  for  all 
i<r . 


5.  Factorization  of  a  finite  band  matrix. 


Theorem  B.  A  TP  nonsingular  (1 ,r )-banded  nxn  matrix  A  can  be  factored  as 


A  =  L 


<*> 


j,(  1  )  D  1  ) 


..  u 


(r) 


with  a  unit-diagonal  (1,0)-banded  TP,  all  k  ,  D  a  diagonal  TP,  and 


a  unit- 


diaqonal  (0,1) -banded  TP,  alJ  k 


Proof.  We  obtain  the  factorization  by  the  standard  device  of  elimination.  In  the 
typical  step,  we  have  a  nonsingular  TP  (i,r)-banded  matrix  A  with  zeros  already  in  band 
-i  in  columns  1,  ...,  k-1  .  From  it,  we  obtain  the  matrix  B  by  subtracting  c  times 


k 

■4 

xxx 
xxx 
0  x  x 
0  x 
Uk  ♦  ft 


x 

X  X 

xxx 

xxx 

xxx 

xxx 


£+k+r 

4- 

x 

X  X 

xxx 
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row  A+k-1  from  row  A+k  .  Thus  B  differs  from  A  only  in  row  A+k  and  there  only  in 
entries  k,  k+1,...,  A+k+r-1  because  of  the  zeros  in  the  other  entries  in  row  A+k-1  .  In 
particular,  B  is  again  (A, r)-banded,  with  zeros  in  band  -A  in  columns  1,...,k-1. 

We  now  choose  c  so  that  also  B(£+k,k)  •  0  .  If  A(£+k,k)  “  0  ,  then  the  choice  c  » 
0  will  do.  Otherwise  A(A+k,k)  y  0  and  then  necessarily  also  A(£+k-1,k)  y  0  (since 
A(A+k-1,k)  »  0  would  throw  its  shadow  to  the  left  because  A  is  nonsingular,  and 
A(£+k,k)  «  0  would  follow).  But  then  the  positive  number 

c  -  A{ A+k ,k)/A( A+k-1 ,k) 

does  the  job. 

Note  that  B  -  CA  ,  with  C  the  identity  matrix  except  for  a  -c  in  position 
(A+k, A+k-1)  .  since  the  action  of  C  is  undone  by  adding  c  times  row  A+k-1  to  row 
A+k  ,  it  follows  that 

A  »  EB 

with  E  the  identity  matrix  except  for  a  c  in  position  (i+k,£+k-1)  . 

In  order  to  carry  out  these  steps  repeatedly,  we  need  to  know 


lam.  B  is  again  TP. 


Proof,  since  B  differs  from  A  only  in  row  A+k  ,  we  only  need  to  consider  minors 
of  B  which  Involve  row  A+k  .  Among  these,  we  only  need  to  consider  those  minors  which  do 
not  involve  row  A+k-1  since  the  others  retain  their  (nonnegative)  value  in  going  from 
A  to  B  .  Thus  we  must  show  that 

b( j )  »  o 

whenever  I  and  J  are  index  sets  of  like  cardinality  and  I  contains  A+k  ,  but  not 
A+k-1  .  Let 

I*  !-  { iel  :  i  <  A+k)  ,  I"  j-  (i£I  s  i  >  A+k}  . 


Then 


cA(I''A+k-1'I") 

>7 
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and  there  is  nothing  to  prove  unless,  as  we  now  assume, 

>  0  # 

which  implies,  by  Hadamard's  inequality,  that  every  principal  minor  of  the  corresponding 
submatrix  is  strictly  positive.  We  must  then  show  that 

>  (-o  • 

For  this,  let,  correspondingly, 

J  -  J'u 

with  J'  the  | I' |  smallest,  and  J"  the  |I"|  largest,  elements  of  J  . 
we  claim  that 


I1  ,i+k-1,I* 


a 


) 


j  ,J 

This  inequality  follows  by  1 1' I -fold  application  of  the  inequality 


-)  • 


c( 

<*>  — 

c( 

valid  for  any  TP  matrix 


1 , 

i,., 

i ,  • 


c  , 


S-1 ,S+1 , . ,t+1 j 
.  .  .  ,  t 

S  ,8+2 ,.,t+1j 
•  «  •  ,  t 


^  f  2 ,  .,s-1 , s+ 1 ,  . , t+ 1  j 
^  v  2 , .  .  .  .,t 

^,2,.,  s  ,s+2,.,t+1\ 
v2, •  •  ,  • ,  t 


because  of  the  identity 


C 1 ,  .,s-1  ,s+1 ,  .,t+1 1  (2,  • ,  s  ,s+2,.,t+1\ 

CM,.  .  .  .,  t  '  c'2,.  .  .  .,  t  ' 

1 ,  • ,  s  ,  s+2 ,  *  , t+ 1  r2,  •  , s- 1 , s+ 1 , . , t+ 1 1 

CM , .  .  .  .,  t  1  *-2,.  .  .  .,  t  J 

[2,...,t+U  r  1,.,s-1,s+2,  .,t+1-| 
t  1  l2,.  .  .  .,  t  1 


valid  for  such  ma+rices.  This  identity  is  proved,  e.g.,  in  Karlin  [7;p.8]  .  It  may  also  be 

proven  by  SDI  applied  to  the  (t+1)x(t+1)  matrix  obtained  by  adjoining  to  the  first  t+1 

T 

rows  and  t  columns  of  C  the  additional  column  (1,0,**«,0)  ,  and  taking 

^  r 1 , • , S~ 1 , 8+2 , • , t+ 1 j 
2,*  .  .  *,t 

as  the  pivot  block. 

Unfortunately,  the  corresponding  argument  involving  dropping  of  the  last  few  rows  and 


columns  reverses  the  sign  in  (*)  and  so  provides  the  irrelevant  inequality 


A(*;k;p/A( 


i+k-1 , 1" 
j 


)  <  a^V^*'1) 
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Instead,  we  observe  next  that 


This  follows  from  the  fact  that 


A(t;k;i;)/A(t+h-1;i:) 


with  the  matrix  C  given  by 


C(y,v)  ;=  a(,^J  ,  all  M,v  , 

V,  J 

hence  TP  (by  SDI),  and  therefore  the  ratio  is  monotone  nondecreasing  in  j  for  j  to  the 
left  of  J”  ,  while  the  strict  positivity  of 

implies,  via  Hadamard's  inequality,  that  a(  +k  1 )  >  0  ,  and  so  k  <  j  (recall  that 

A(i+k_1  1-0  for  v  <  k  ) . 

V  V  J 

This  leaves  us,  finally,  with  the  task  of  showing  that 

But  that  is  now  obvious  since  A(i,k)  «  0  for  i  >  4+k  ,  hence 


A(4;k)A(fj 


A(4;ki/A(4+k-1) 


in 


we  conclude  that  a  nonsingular  ( 4 ,r)-banded  TP  matrix  A  can  be  factored  as 


A  -  B 


( k) 

with  B  again  TP  but  only  (4-1 ,r)-banded  ,  while,  for  each  k  ,  E  is  the  identity 


matrix  except  for  some  nonnegative  c^  in  position  (4+k,4+k-1)  •  But  then 


(4) 


,o>. 


.( n-4 ) 


is  a  ( 1 ,0 )-diagonal  matrix  with  unit  diagonal  and  the  nonnegative  number  c^  in  position 
(4+k,4+k-1)  ,  k»1,...,n-4  ,  and  zero  everywhere  else.  Consequently,  is  (1,0)-banded 

and  TP. 


we  conclude  that  a  nonsingular  ( 4, r) -banded  TP  matrix  A  can  be  factored  as 


r<4) 


..L(,)B 


-1 


1  - 


t 


with  B  a  (0,r)-banded  TP  matrix  and  each  a  TP  (1,0)-banded  matrix  with  unit 

diagonal.  Applying  this  last  statement  to  BT  and  transposing  the  result  finishes  the 
proof  of  Theorem  A.  I | | 

6.  Factorization  of  a  (bi) infinite  band  matrix. 


Theorem  C.  A  TP  (bi) infinite  m-banded  matrix  A  whose  rows  and  columns  are  linearly 
independent  can  be  factored  as 

A  * 

/  ki 

with  each  r'  a  TP  one-banded  matrix  with  maximum  entry  in  each  column  equal  to  1  and 
D  a  TP  diagonal  matrix  with  0  <  D(j,j)  <  max^  A(i,j),  all  j  . 


Proof.  If  A  is  biinfinite,  then  we  know  from  Theorem  A  that  all  principal  sections 
for  some  band  of  A  are  nonsingular.  Assume  without  loss  of  generality  that  the  zeroth 
band  is  such  a  distinguished  band  and  that  A  is,  more  explicitly,  ( l ,r )-banded.  Then  we 
know  that,  for  all  n. 


Ajj  i*  A(-n,  n] 

N*N 

is  nonsingular.  If  A  is  only  infinite,  A  £  R  say,  then  we  know  from  the  Corollary  to 
the  Proposition  in  Section  4  that,  for  all  n, 

An  s  a  A(  1 , . . .  ,n) 


is  nonsingular. 

In  either  case.  Theorem  B  assures  us  that  A„ 


.  -  i(1)  r(1)  n  „(1> 

A  —  L  ...X,  DU 
n  n  n  n  n 


has  a  factorization 
(r) 


U 


with  i/*^  unit-diagonal  (1,0)-banded  TP,  unit-diagonal  (0,1)-banded  TP,  and  D 

n  n  n 


diagonal  TP.  We  intend  to  let  n  go  to  infinity  and  therefore  must  deal  with  the  fact  that 


these  factors  may  not  be  bounded  independently  of  n 


For  this,  define  one-banded  matrices  S(  ** ,  ....  S(r*  as  follows.  Starting  with 
M  1  ,  define  S  ,  ...,  S  successively  by 

with  the  diagonal  matrix  having  max^M*  )  (  i,  j  )  in  its  j-th  diagonal 

position.  This  number  cannot  be  zero  {by  induction  on  X  )  since  M*  * 1  »  1  and  l/k3 

n 

unit-diagonal.  It  follows  that  each  S(_,t*  is  a  (1,0)-banded  TP  matrix  with  maximum 


is 


entry  1  in  each  column,  and 

,  (*> 


t<1>  e(-i) 

L  •  •  •  L  *  S 

n  n 


s'-V0* 


Now  continue  the  process,  starting  with  M* 1  *  s»  M*°*Dn  »  getting  successively  S*1*,  ..., 
S(r)  by 

S(X,  s,  M(k,u(ic)(jl<k+1,)-1  / 

with  M^+1>  the  diagonal  matrix  having  max^M^u' )  ( i,  j )  in  its  j-th  position. 


We  arrive  at  the  factorisation 

.(-*>  S<_1,Sn)  S(r)Mtr)  -  R(1) 

e  e  e  O  O  t  *  *  O  H  •  R  •  e  * 


s 


(m) 

r'  ;E 
n  n 


,00 


with  each  one-banded  TP  and  maximum  entry  1  in  each  column,  and  En  a  diagonal  TP 

matrix.  We  claim  that 


0  <  E  <j,j>  <  max  A(i,j)  ,  all  j 

"  i 


We  know  that 


A(i,j) 


r  ...  J  R^’u.j.)  Rl2)(j.«j,)  .. 

.  n  i  n  i  4 

31  3m 


C'3 


with  all  summands  nonnegative.  Further,  for  at  least  one  choice  of  i  ,  one  of  the  summands 
is  just  En(j,j)  since,  starting  with  jB  “  j  ,  we  can  pick  3m-2’  •••»  “!  *  in 

sequence  so  that  < l*-! '3**  “  1  ’  But  then  *  En<3,3>  * 

We  can  now  let  n  go  to  infinity  through  a  subsequence  of  ■  in  such  a  way  that  each 

necessarily 


of  the  matrices  R*k'  converges  entrywise  to  some  (bi)lnfinite  matrix  R* * ^ 
n 


one-banded  TP  with  maximum  entry  1  in  each  column,  and  En  likewise  converges  to  some 
diagonal  matrix  0  satisfying  0  <  D(j,j)  <  max^A(i,j)  ,  all  j  ,  while 


R<m>  D 
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But  then  0  <  D(j,j)  .  all  j  ,  since  otherwise  A(*,j)  »  0  ,  contradicting  the  linear 


independence  of  the  columns  of  A  .  Ill 


Corollary.  Let  A  be  a  (bi) infinite  TP  m-banded  matrix  whose  rows  and  columns  are 
linearly  independent.  Then  A  is  the  limit  of  strictly  m-banded  (bl) infinite  TP  matrices 
and  thi3  limit  is  uniforr  (i.e.,  in  norm)  if  A  is  bounded. 


Proof.  Replace  each  zero  entry  in  the  two  interesting  bands  of  r'k'  above  by 

( to 

e  >  0  to  obtain  the  strictly  one-banded  TP  matrix  R  ,  all  k  .  Then 

A  R(1)  ...  R(m>  D 

c  e  e 

is  strictly  m-banded  TP  (as  a  product  of  strictly  banded  TP  matrices)  and  converges 
entrywise  to  A  as  e  — *  0  .  Since  the  entries  of  R*’5^  are  bounded  by  1  while  those 


of  D  are  bounded  by  lAl^  ,  this  convergence  is  obviously  uniform  in  case  lAl^  <  “  .III 


Remark.  The  assumption  that  the  rows  and  columns  of  A  are  linearly  independent  is 
not  bothersome  in  the  intended  use  of  this  corollary  in  fl)  since  there  A  is  even 
boundedly  invertible.  But  it  would  be  nice  to  know  whether  this  assumption  is  necessary.  We 
note  that  Metelmann  [8]  has  obtained  strictly  one-banded  TP  factorizations  for  finite 
strictly  (f.r)-banded  TP  matrices  and  that  Cryer  [4]  has  obtained  one-banded  TP 
factorizations  for  arbitrary  finite  TP  matrices.  But,  the  procedure  given  by  Cryer  may 
produce  more  than  m  1-banded  factors  unless  the  matrix  is  strictly  m-banded. 
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